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Abstract. We formulate a general and exact method of partial wave decomposition (PWD) of any nucleon- 
nucleon (NN) potential and any three-nucleon (3N) force. The approach allows one to efficiently use 
symbolic algebra software to generate the interaction dependent part of the program code calculating the 
interaction. We demonstrate the feasibility of this approach for the one-boson exchange BonnB potential, 
a recent nucleon-nucleon chiral force and the chiral two-pion-exchange three-nucleon force. In all cases 
very good agreement between the new and the traditional PWD is found. The automated PWD offered 
by the new approach is of the utmost importance in view of future applications of numerous chiral N3LO 
contributions to the 3N force in three nucleon calculations. 

PACS. 21.45.-v few-body systems - 21.30.-x nuclear forces - 21. 45. Be two-nucleon system - 21.45.Ff 
three-nucleon forces 



1 Introduction 

The standard way to set up calculations of two- and three- 
\ nucleon systems is a partial wave decomposition (PWD). 
Especially at low energies, i.e. below the pion production 
threshold, this procedure is still most commonly used, 
despite the advent of the approaches which use a direct 
three-dimensional notation [1,2,3,4]. 

Recently we proposed a formulation of the two- and 
three-nucleon system 5,6 , which is based on scalar spin- 
momentum operators and accompanying scalar functions 
depending only on the momenta of the system. This for- 
mulation is based on the most general operator structure 
a nuclear force given in momentum space can have. The 
two- and three-nucleon equations are obtained by carry- 
ing out traces over the spin-momentum operators building 
the nuclear force. The same approach can be used to ob- 
tain partial wave projected matrix elements of the poten- 
tial and the transition operators. Taking traces of spin- 
momentum operators lends itself to the use of symbolic 
algebra software to obtain general expressions for those 
matrix elements. 

In this work we demonstrate that by algebraic oper- 
ations general expressions for the partial wave decompo- 



sition of any nucleon-nucleon (NN) potential can be ob- 
tained. This method will be presented in Sec. [H and ex- 
plicit expressions for calculating specific matrix elements 
of NN potentials are given in an appendix. Numerical com- 
parisons between our suggested methods and the stan- 
dards partial wave decomposition using recent chiral NN 
forces [7] and the one-boson-exchange NN potential BonnB 
[8] are presented in Sec. [3l In Sec. [4] we demonstrate that 
the same method can be extended to the PWD of a three- 
nucleon force. As a numerical example given in Sec. [5] we 
take the two-pion-exchange (TPE) chiral NNLO three- 
nucleon force. Finally we conclude in Sec. [6] 

2 Partial wave decomposition of the NN 
potential 

We start, as in Ref. [6], by projecting the NN potential 
on the two-nucleon (2N) isospin states \ {W)tm t ) = \tmt). 
Furthermore, we assume that there is no isospin mixing 
but allow for charge independence and charge symmetry 
breaking, and thus for a dependence on to*: 

(t'm' t \V\tm t ) = 5 t , t 8 m . tmt V tm \ (1) 
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It is well known that the most general form of the NN 
force, which is invariant under rotations, parity and time 
reversal can be expressed by six scalar spin-momentum 
operators [10], which we choose as 

u>i(o-(l),o-(2),p',p) = 1 , 

w 2 (<T(l),cr(2),p',p) = a{l)-a{2) , 

w 3 (ct(1), cr{2),p',p) = i (rr(l) + rr(2)) • (p x p') , 

» 4 Wl),a(2),))» = a(l).(pxp>(2).(pxp') , 

w 5 (a(l) ) <T(2),p\p) = <T(l)-(p'+p)<r(2)-(p'+p) , 

w 6 (cr(l),a(2),p',p) = tj{l) ■ {p' - p) a{2) ■ {p' - p) . 

(2) 

Thus the isospin projected potential can be expresses as 

6 

V tmt (p',p) =]r^(p',pH(<x(l),<x(2),p',p). (3) 

The expansion coefficient here are scalar functions fi(p', p) 
that depend on two vector momenta p' and p, more specif- 
ically on the magnitudes of the vectors and and the cosine 
of the relative angle between them. In order to determine 
the functions ft, we evaluate the spin dependence analyt- 
ically by taking traces with the operators of Eq. ([2]) and 
thus arrive at a system of six coupled linear equations 

6 

Y Tr (wi Wj ) ft = Tr (V tmt Wl ) , i = 1, 2, . . . , 6, (4) 
i=i 

which has a unique solution provided that p' ^ p. (When 
p' = p only five out of the six operators u>i are sufficient, 
since u>2 is linearly dependent on W4, W5, and Wq-) 

Our task is to obtain matrix elements of V trrlt in the 
basis of states | p(ls)jrrij), where the relative angular mo- 
mentum I and the total spin s are coupled to the total 
angular momentum j with its projection rrij. When calcu- 
lating NN observables one usually sums angular momenta 
j up to a certain j max at which the calculation is con- 
verged. For calculations of NN observables below the pion 
production threshold convergence is reached for j max < 8. 

To obtain the potential matrix element in the basis 
I p(ls)jrrij), the four-fold integral 

(p'(l' S )j mj I V tm < I p(ls)j mj ) = 

x ^2c(l,s,j;mi,mj - mi,mj) Y^^O', <f>') Y lmi {6,(j>) 

mi 

x(srrij — m'i \ V tmt (p',p) \ srrij — mi) (5) 

needs to be evaluated. Here c(Z, s, j;mi,rrij — mi,rrij) are 
the standard Clebsch-Gordan coefficients, and Yi mi (9,<f) 
the spherical harmonics calculated for the angles corre- 
sponding to the directions of the momenta p' and p. This 
quantity does not actually depend on rrij, so instead of 
Eq. ([5]) we can calculate 



H(l',l,s,j) = - - (P'( l ' s )i m j I V I p{ls)j mj ). 

(6) 

Since now the integrand is a scalar, it is possible to reduce 
the number of integrals to one: 

r 1 

H{l',l,s,j) = ^ 2 / d{cose')j-^ 

_1 rrij=-j 

I' 

x Y c ( l ', s ,3'-, m 'i, m 3 - m 'l, m j) 

m> = -l> 
I 

x Y / c(l,8,j;mi,m,j - mi,rrij) 

mi=—l 

x Y Vm , (^',0)^,(0,0) 

x(srrij — m'i \ V(p',p) \ srrij — mi). (7) 

In the NN system we can choose the z-axis to be 

p=(0,0,p). (8) 

Then the direction of the vector p 1 is given as 

p' = O'sin0',O,p'cos0'). (9) 

Most importantly however, the matrix element in the 2N 
spin space, 

(s m.j — m[ \ V(p' ', p) \ s mj — mi) = 

6 

(smj -m' l I YMp'>P' x ) 

i=l 

w i (a{l),cT(2),p , ,p)\sm j -mi) (10) 

can be calculated analytically Q- The three sums in Eq. ([7]) 
over mj, m\ and mi can be written out explicitly so the 
integrand over x = cos 9' can be prepared once for all 
NN potentials that are represented in the form of Eq. © . 
The resulting expression is given in terms of p, p' , x and 
the expansion coefficients f(i) = fi(p',p,x). As example, 
evaluating Eq. ([7]) for the 1 5'o channel leads to 

1 

H(0, 0, 0, 0) = 2tt j dx (/(l) - 3/(2) + f(4)p 2 p' 2 (x 2 - 1) 
-1 

~f(5)(p 2 +p' 2 + 2pp'x) 

-f(G)(p 2 +p' 2 -2pp'x)). (11) 

Further examples are collected in Appendix lAl 



For this calculation symbolic software like Mathematica 
© [II] proves very useful. Here it is particularly simple with 
the concept of the Kronecker product which makes such matrix 
elements simple matrix elements in the four-dimensional space. 
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3 Application of the new method to the 2N 
potential 



In this section we want to give two examples a PWD of an 
NN potential based on our new method and compare with 
results obtained in the traditional way. First we consider 
the one-boson exchange potential BonnB [8] . Here we use 
the operator form of this potential presented in [6]. In 



Figs. [TH6] we show matrix elements (p'(l's)jr 



ytm t 



Pa(ls)jm,j) for several partial waves for a fixed value of 
Po = 1 fm _1 as a function of the momentum p' . In Figs. 1 
through 3 selected partial wave projected potential ampli- 
tudes for the BonnB potential are shown. The agreement 
of the matrix elements calculated with our operator based 
formulation and the original PWD based on a helicity for- 
mulation [9] is excellent. That is not too surprising, in 
both cases the numerical accuracy is determined by the 
integration over cos6>'. 

As second example we take the neutron-proton version 
of the chiral NNLO potential [7] and use the same set of 
the parameters as given in [6]. In Figs. 0H6] selected partial 
wave projected potential amplitudes are shown for this 
potential. Again, the agreement between the two methods 
is excellent. 



4 Extension to Three- Nucleon Forces 



The arguments used in Sec. [2] for the 2N potential can be 
applied to any three-nucleon (3N) force. Now we use the 
3N states | pqf3) [12] in the so-called LS-coupling 



pq/3) =\ pq(l\)L(s±)S(LS)JMj) \ (t\)Tm T ) , (12) 



where the quantum numbers for the relative angular mo- 
menta I (within the pair (23)) and A (between the pair (23) 
and nucleon 1) are coupled to the total angular momen- 
tum L. In the spin space the spin of the (23) pair is coupled 
with the spin \ of nucleon 1 to the total spin S. Finally L 
and S are coupled to the total 3N angular momentum J 
with the projection Mj. The isospin 3N state, where we 
couple the total isospin of the (23) subsystem t with the 
isospin i of the third nucleon to the total 3N isospin T 
with the projection mr, is constructed in the same way 
as the 3N spin state. Since we want the states | pq(3) to be 
antisymmetric with respect to the exchange of particles 2 



l+s+t 



and 3, we require additionally that (—1) 

We start with the eight-dimensional integral 



-1. 



{p' q \l l \')L\s ll -)S'{L'S')JM J \ 



V 3N | pq(l\)L(s-)S(LS)JMj) 



dp 'J dq 'J dp j dq 



c(L', S', J; m L ,,Mj - m L ,,M.,) 

m L , 

53 C ( L > S > J \ m L, M.j - m L ,Mj) 

(( s l)S'Mj-m L , | V 3N (p',q',p,q) \ ( s ±)SMj-m L ), 

(13) 

where 

i 

yt'™ L (p> = 53 c ^> \L; m h m L - mum L ) 

m L = -l 

(?)• (14) 

This quantity is independent of Mj, and instead of Eq. |T 
we can calculate 

G(l' , A', L , s , S , I, A, i, s, S, J) 
^^tVi £ (p'q'(l'X)L'(s'^)S'(L'S')JMj | 



Mj=-J 



V 3N | Pq (lX)L( s ^)S(LS)JMj) 



(15) 

Since the integrand is now a scalar, it is possible to re- 
duce the number of integrals from eight to five. Namely 
we take first p = z and then consider all scalar products 
among the p, q, p' and q' unit vectors. They depend on 
the following quantities: 8 q , p i , 8 q i , <p p > — <p q , <j) q i — <t>„ and 



4>p' ~ 4>q' ■ Since the last three are not independent, it is 
possible to set additionally <f) q = 0. Of course, taking p = z 
and <j) q — is only one possibility. The best choice (from 
the computational point of view) of the five integration 
variables might depend on the form of V 3N (p' ,q' ,p, q), 
which is an operator in the 3N spin space. This quantity 
G has to be multiplied by the matrix element of isospin 
operator I 



{{t' l -)T'm TI \i\{t l -)Tm T ), 



(16) 



which can be worked out exactly, independent from the 
momentum and spin spaces. 



5 Numerical example for a 3N force 

We consider just one example of the 3N force: a two-pion- 
exchange contribution to the chiral NNLO 3N force as 
given in [T3] , V^p^ ■ In particular we take the part of V^p^ 
which is symmetric under the exchange of nucleons 2 and 
3 
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Fig. 1. Comparison of the traditional (crosses) and the new (solid line) method of PWD for the Bonn B potential. The matrix 



elements of selected uncoupled channels with t 
as a function of the p' momentum. 



1 ( l S , 1 D 2 ) and t = ( 3 G 4 , 3 H 5 ) are shown for a fixed value of po= 1 fm 



V {1) = F\Oi ■ q 2 cr 3 • q 3 r 2 ■ r 3 

+ F 2 cr 1 ■ (q 2 x q 3 ) cr 2 q 2 (r 3 - q 3 n • (r 2 x r 3 ) , 

(17) 

where q i = p' ' i — p i and p i (p'J is the initial (final) mo- 
mentum of nucleon i. Further, gA, and F„ refer to 
the nucleon axial vector coupling constant, pion mass and 
decay constants, respectively, while Ci are low-energy con- 
stants from the subleading pion-nucleon Lagrangian. The 
scalar functions F\ and F 2 are 



Fx = 



9A 



2i*7r / (q 2 2 + 77l2)(g 3 2 +m 2) 



4cim 2 2c 3 
1 ~ p^+p2^-Qs 



and 



Fa 



.9-1 



1 



<~4 



(18) 



(19) 



.2F w y ( < 7 2 2 +m2)(q 3 2 +m 2 ) F 2 

We consider first the two isospin matrix elements 

Ix{t',T',mT>,t,T, m T ) 

= ({t'\)T'm T ,\T 2 -T 3 \{t~)Tm T ) (20) 



and 



I 2 {t 1 , T', m T > , t, T, m T ) 

= {{t'\)T'm T , | n ■ (r 2 x t 8 ) | {t\)Tm T ). (21) 



The matrix elements Ji(t' , J 1 ' , mr' ,t,T,m,T) are particu- 
larly simple and are given as 

h{t',T',m T ' ,t,T, m T ) 

= (2t(t + 1) - 3) St,t> S mT , mT , ■ (22) 

The (purely imaginary) matrix elements 
l 2 (il ', T', mT' j T, %) can be written as 

I 2 (t',T',mT>,t,T, m T ) 

= iVs (-l) t+1 St + t',i S T i S T , i S mTt m T ,. (23) 

From Eqs. (f2"2")l and (|2"3"|) we defer immediately that the two 
parts of the considered 3N force (fT7|) will not contribute 
simultaneously to the same matrix element 
G(/', A', L 1 , s', 5', A, -Z^, s, 5*, */)■ 

For this first simple study we construct 16 | (3 ) states 
for J = ~ and positive parity 7r = (— 1)' +A satisfying the 
additional condition I < 2 and A < 2. Their quantum 
numbers are given in Table [T] 

Next we perform the steps described in Sec. 0] and 
obtain 256 integrands G(l', A', L', s' f S', I, A, L, s, S, J) = 
G(J3',p) such that 

G((3',(3) = ({t' l -)T' mT , | (p'q'p I V m I pq(3) | {t l -)Tm T ) 

(24) 

Here we show just few (relatively simple) examples. 
(Note that in the first two cases the scalar nature of the 



dp' dq dp dqG(l3',(3) . 
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Table 1. List of | /?) states for J n = ± , / < 2 and A < 2. 
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Table 2. Stability of the selected five fold integrals with re- 
spect to the number of Gaussian points N (see text). In the 
last row the values obtained with standard PWD [14II15| are 
given. All matrix elements are in fm . 



N 


G(l,l) 


G(2,l) 


G(6, 12) 


G(5, 11) 


12 


443.565 


1200.160 


-5.52616 


-5.24720 


24 


443.618 


1200.223 


-5.49311 


-5.48527 


36 


443.618 


1200.219 


-5.49290 


-5.48630 


48 


443.618 


1200.219 


-5.49290 


-5.48626 


standard 


443.618 


1200.219 


-5.49274 


-5.48597 



PWD 



Table 3. Stability of the selected six fold integrals with respect 
to the number of Gaussian points N (see text). All matrix 
elements are in fm J . 



G(/3',/3) functions is clearly visible.) 
1 



G(l,l) 
G(2,l) 



16tt 2 

i 



F 1 h(l,T\m T/ ,l,T,mT)q 2 - q 3 , 
F 2 l2(0,T',mT>,l,T,m T ) 



167T 2 %/3 

1 



G(5,ll) = ^=F 2 I 2 (0,T , ,m T ',l,T,m T ) (q 2 -q 3 ) 

x y 2 , 2 (Pi q) (nx - iqiy) yll'^ip', q') 

+2q 4z yl 1 { °(p',q') 
-V2(q ix +iq4y)yti 1 (p 1 q)j , 
G(6,12) = \F x i x {l,T' ,m T ,,l,T,m T )y^\p' A 1 ) 

2(fe(<73x + m y ) 
~(q2x + iq2 V )q3z)y 2 ' i2 1 (p, q) 

+ 2«((?2yfe - q2xq3y)y 2 ' 2 (P) <?) 



x iq3y) 

-(q2x - iq2 V )q3z)y 2,2 (p,q)) , 



(25) 



where q 4 = q 2 x q 3 and the Cartesian components of 
the q i vectors are denoted as g, x , qi y and qi z . Of course, 
vectors q t are now expressed in terms of the initial and 
final Jacobi momenta in the following way 



q -q 
, i 



<? 2 = P - 2« 



P~- 2 q 



13 



P 



P - 2 q 



(26) 



Many of the G(/3', (3) functions are very lengthy. Using 
symbolic manipulations software like Mathematica © [llj 



N 


G(l,l) 


G(2,l) 


G(6, 12) 


G(5,ll) 


12 


443.510 


1200.365 


-5.55695 


-5.14678 


24 


443.619 


1200.218 


-5.49277 


-5.48785 


36 


443.618 


1200.219 


-5.49290 


-5.48626 


48 


443.618 


1200.219 


-5.49290 


-5.48626 



one can write down these expressions directly as a part of 
a Fortran or a C code. The expressions given above have to 
be integrated over five angles. With the choice of variables 
discussed in Sec. [4] we calculated the five fold integrals 
for fixed magnitudes of the momenta p, q, p' and q' . They 
were chosen quite arbitrarily: p— 1 fm , q— 2 fm , p'= 
3 fm -1 , q'= 4 fm -1 . For the isospin quantum numbers we 
assumed that T = T' = tut = tut 1 = \- The 3N force 
parameters were taken as Qa= 1.29, F n — 92.4 MeV, 
138.0 MeV/c 2 , c x = -0.81 GeV" 1 , c 3 = -3.4 GeV" 1 , c 4 = 
3.4 GeV -1 . Such single five fold integrals can be easily 
calculated on a PC. For a simple test we chose the same 
number of Gaussian integral points N in each of the five 
dimensions and checked the stability of the integrals with 
respect to N . Even in this most simple way, we found a 
good convergence, which is shown in Table [2l To check nu- 
merically the scalar nature of our integrands we calculated 
also the corresponding six fold integrals, where (j> q was an 
additional integration variable and checked their stability 
(see Table [3]). We see from Tables [5] and |3] that the results 
from the five and six fold integrations agree very well with 
each other. We calculated the same matrix elements using 
standard PWD [14 , 15j . The results are shown in the last 
row of Table [2l They agree very well with the numbers 
obtained using the new method. Thus we conclude that 
the presented procedure might be used to efficiently auto- 
mate the very cumbersome standard PWD of the 3N force. 
Calculation of all matrix elements for the whole grids of 
p, q, p' and q' points, especially for bigger values of J, will 
require an implementation of the algorithm on parallel 
architectures. Implementing such a code is much simpler 
than implementing a code based on the traditional PWD. 
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6 Summary and Outlook \ ( 

tf(l,l,l,0)=27r J dxl [/(l) + /(2)] 

We propose a new general method of calculating the PWD -i 

of any NN potential which is given in momentum space +2/(3) pp'(x 2 — 1) + /(4) p 2 p' 2 x(l 

in terms of operators acting in the 2N spin and isospin 

spaces. Once the expansion coefficients in the operator ~ [2pp + x(p + p )J 

basis are determined, the matrix elements of interest are / 1 

obtained as simple one dimensional integrals, independent +/(6) [2pp — x(p +p )J > . 

from the particular form of the NN potential. We demon- ' 
strate this method with two examples, the one-boson- 
exchange BonnB potential and a recently derived chiral f unc tions for j = 1 are given by 
NNLO potential. In both cases the agreement between our 
new method and traditionally employed one to obtaining 

partial wave projected matrix elements is perfect. i r 

This method can be extended to treat three-nucleon HiX 1 1) = 2tx \ dx < [/(l) + 3/(2)] 

forces, for which a traditional PWD is a formidable task. J | 

In order to show the power of our new method, we per- 
formed a simple feasibility study for the TPE chiral 3N 



-/(4) p 2 p' 2 x(x 2 - 1) 



force. A very good agreement between the new method — /(5) % \2p'px + (p' 2 + p 2 )] 

and standard PWD shows that the very complex tradi- 
tional approach can be replaced by this straightforward +/(6) x [2p'px — (p 12 + p 2 )] 
and efficient algebraic method. An 'automated' PWD for 

the 3NF is of great importance in view of awaiting ap- 1 
plications in the 3N systems. The numerous contributions 
to the chiral N3LO 3N force pTBlTT?] can be handled only 

with such an efficient and simple procedure. -l 

+/(3) p'p(x 2 - 1) 



#(1,1, 1,1) = 2^ / dxl [/(l) + /(2)]x 
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tly performed on the supercomputer cluster of the JSC, 
Julich, Germany. 



A Integrals for the 2N potential 



+/(6) [(3a; 2 - l)p' 2 + 2p 2 - Ap'px] | 
l , 

#(2,0, 1,1) = 2^ J dx\f{A)p' 2 p 2 { 



This appendix collects the expressions for the partial wave , . > > , 

decomposition of the NN potential for the partial waves [2p + [x — l)p + 4p px\ 

with total 2N angular momentum j < 2. The expressions 

for higher j can be obtained on request from one of the + f(6) \2 12 + (x 2 — 1) 2 — 4 ' x] 

authors (JG) either in analytical form or as Fortran or C w L 7 j 

code. 

The second function for j — is {H(Q, 0,0,0) was given 

in Sec. [2|): and 
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#(2,2,1,1)= f/<fcj 



3[/(l) + /(2)] (3a; 2 -1) 



+18/(3) p'px(x 2 -l) 
+/(4) p'V(14x 2 -5-9a; 4 ) 
+/(5) [(l-3x 2 )(p' 2 + p 2 )-Ap'px\ 

+/(6) [(1 - 3a; 2 )(p' 2 + p 2 ) + VH | • (33) 
The functions for j = 2 are given by 

H(2,2,0,2)=7r J dx(l-3a; 2 )j [-/(l) + 3/(2)] 

-/(4)p'V(^ 2 -l) 
+/(5) [(p' 2 +p 2 ) + 2p'px] 

+/(6) [(p ,2 +p 2 )-2p'px] 1 , (34) 



Ji(3, 1,1,2) = 2tt^ | | /(4) p' 2 p 2 ^ 2 - 1) 
+/(5) [2p' 2 +p 2 x(5x 2 - 3) + 2p'p(3a; 2 - 1)] 

+/(6) [2p' 2 +p 2 x(5x 2 - 3) - 2p'p(3a; 2 - 1)] j , 

(38) 



if (3, 3,1,2) = - y da; |5a;(5x 2 - 3) [/(l) + /(2)] 

+10/(3) p'p( 1 + 5a; 4 6a; 2 ) 

+/(4) p' 2 p 2 a;(44a; 2 - 25a; 4 - 19) 

+/(5) [(p' 2 +P 2 )x(3 - 5x 2 ) + 2p'p(l - 3a; 2 )] 



+/(6) [{p 12 + p 2 )x(3 - 5a; 2 ) - 2p'p(l - 3a; 2 )] 



(39) 



if (2, 2,1,2)=* y dx| [/(l) + /(2)] (3a; 2 - 1) 

+2/(3) p'px(x 2 - 1) 

+/(4) p'V [^ 4 ~ 2a; 2 + l] 

+/(5) [(3a; 2 -l)(p ,2 +p 2 ) + 4p'pa; 3 ] 



+/(6) [a;(3a; 2 - l)(p' 2 + p 2 ) - Ap'px A \ j , 

ff(l,l,l,2) = y y da^x [/(l) + /(2)] 

+5/(3) j/p(l - a; 2 ) 

+2/(4) p'V^(l - a; 2 ) 

+/(5) [a;(p' 2 +p 2 )+p'p(3a; 2 -l)] 

+/(6) [a;(p' 2 +p 2 )-p'p(3a; 2 -l)] | , 
l , 

if (1, 3, 1, 2) = 2tt^ J dx J /(4) p' 2 p 2 ^ 2 - 1) 
+/(5) [p' 2 a;(5a; 2 - 3) + 2p 2 a; + 2p'p(3a; 2 - 1)] 

+/(6) [p /2 a;(5a; 2 - 3) + 2p 2 a; - 2p'p(3a; 2 - 1)] 



(35) 



(36) 



(37) 
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Fig. 3. The same as in Fig. [2] for one selected coupled channel case with t = 1. 



J. Golak et al.: A new way to perform partial wave decompositions of few-nucleon forces 



9 



1 s 




o 





-0.0005 - 


CM 

E 


-0.001 - 


-0.0015 - 


A 
o 

CL 
> 


-0.002 - 
-0.0025 - 
-0.003 - 


"a. 

V 


-0.0035 - 




-0.004 - 



-0.0045 1 1 1 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 4 4.5 

P' [fm _1 ] 

Fig. 4. The same as in Fig. [1] for the example 
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a chiral NNLO potential (see text for details). 
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Fig. 5. The same as in Fig. [2] for the selected chiral NNLO potential. 
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Fig. 6. The same as in Fig. [3] for the selected chiral NNLO potential. 



